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PLANAR PARTICLE/DROPLET SIZE 
MEASUREMENT TECHNIQUE USING 
DIGITAL PARTICLE IMAGE 
VELOCIMETRY IMAGE DATA 

This invention was made with government support under 
Grant No. NGT 3-52352 NASA Glenn Research Center and 
Grant No. NAG3-2110 NASA Glenn Research Center. 

BACKGROUND OF THE INVENTION 

The present invention relates to particle sizing. It finds 
particular application in conjunction with two-phase flows 
and will be described with particular reference thereto. It 
will be appreciated, however, that the invention is also 
amenable to other like applications. 

Two-phase flows are used in many fuel combustion 
processes such as those found in gas turbine combustors, 
coal furnaces for power generation, and diesel engines. The 
efficiency of two-phase flows directly impacts the efficiency 
of the associated combustion process. For many industrial 
and fossil fuel energy processes, both the droplet/particle 
size and spatial distribution are of interest. The useful 
diagnostic tools should be able to make in situ measure- 
ments without disturbing the flow field so that the measure- 
ments can be meaningful for understanding these flows. This 
suggests the need for a novel non-intrusive optical technique 
that can provide instantaneous measurements of particle size 
and velocity at multiple spatial points in planar (2-D) fields 
so that estimates of the mass flow are obtained. 

Digital Particle Image Velocimetry (DPIV) is a technique 
for obtaining planar measurements of particulate seeded 
flow fields. Light from a pulsed laser is formed into a thin 
sheet to illuminate a planar cross section of a flow. A CCD 
camera is used to record the light scattered by the particles 
at the two (2) instants that the light sheet is pulsed. The fluid 
velocity is determined by analyzing the recorded particle 
image data. Typically, cross-correlation data analysis is used 
to reduce the recorded image data to determine the fluid 
velocities. In most instances, the flow field is artificially 
seeded with tracer particles in order to measure the flow 
velocity. In two-phase flows, the second phase material 
provides the scattering sites for estimation of the flow 
velocity, or at least the velocity of the second phase. The 
objective of this work is to determine the feasibility of 
estimating particle size from Particle Image Velocimetry 
(PIV) image data. Particle size information can most likely 
be extracted provided the imaged particle size exceeds the 
optical system blur circle. The other major factors affecting 
the accuracy to which particle size information can be 
extracted are the optical system f/number, image system 
pixel resolution and dynamic range, the optical properties of 
the particles, and the characteristics of the scattered light 
from the particles. 

Many different optical techniques for making in situ 
measurements of particle/droplet sizes are known. 

An in-line holography system has been used to record 
holographic images of the light scattered from particles in a 
fluid. The hologram has been shown to be that of a screen 
with a circular aperture representing the particle image. 
Holograms have been reconstructed by passing a laser beam 
through the hologram and recording the reconstructed image 
at an on-axis observation plane. The reconstructed image 
contains the Fraunhofer pattern of the imaged particles. The 
size and shape of the particle is then determined by direct 
observation of the reconstructed image. Holography has 
been easily extended to velocity measurement by using 
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double exposure holograms and analyzing the double aper- 
ture type of fringe structure in the reconstructed image to 
calculate the particle separation. 

Laser interferometry has been used to characterize par- 
5 tides using the spatial frequency of the far-field fringe 
pattern in the forward-scatter region. Ovryn 4 uses partially 
coherent light and forward scattering to obtain Poisson Spot 
images of particles suspended in a solution to determine 
their 3-D velocities. Although not a sizing technique, this 
10 work illustrates the information content in the diffraction 
rings surrounding coherently illuminated particles and such 
rings are used to determine additional properties about the 
particle. The difficulty with these holographic and interfero- 
metric techniques lies in the requirement of a high spatial 
resolution detection media, typically photographic film. 

Phase Doppler Particle Analysis (PDPA) is an existing 
technique that is able to make point wise velocity and size 
measurements simultaneously. Droplet sizes are obtained 
using PDPA, which relates the droplet size to the phase shift 
20 of light refracted through the drop and scattered to different 
positions on the receiving lens. The technique known as 
Laser Doppler Anemometry (LDA), which is a technique for 
obtaining point wise velocity measurements, has been modi- 
fied for applications to particle sizing as well. The particular 
25 modifications relate the particle diameter to the LDA signal 
visibility. Other variations of the technique relate either the 
scattered intensity or the phase shift of the scattered light to 
the particle diameter. 

In diffraction-based sizing techniques the particle inter- 
action with the illumination is assumed to be analogous to 
the interaction of the same illumination with a uniformly 
illuminated circular aperture. The inherent problem in dif- 
fraction methods lies in the assumption that a particle acts as 
a uniformly illuminated aperture, which ignores many other 
scattering effects present in light scattering from a spherical 
particle, such as specular reflection, interference and refrac- 
tion effects. Therefore, the Mie scattered signal has been 
used for determining particle characteristics. 

40 The technique of Laser Sheet Dropsizing (LSD) has been 
used to determine droplet diameters. LSD uses the Mie 
signal, which is proportional to the particle diameter 
squared, and the Laser Induced Fluorescence (LIF) signal, 
which is proportional to the particle diameter cubed, for 
45 determining particle size. 

Mie scattering has been investigated from small particles 
in the regime where the particle image is dominated by 
diffraction to determine the minimum particle diameter that 
is detectable for use in pulsed laser velocimetry techniques. 
50 Taking the Fourier transform of the Mie scattered electric 
fields from spherical particles results in multiple peaks in the 
intensity profile, the characteristics of which can be related 
to the particle size. 

The ratio of the projected area of the reflection glare spot 
55 to the cross-sectional area of the droplet has been found to 
be proportional to the square of the ratio of the aperture 
radius to the distance from the droplet to the sensor surface. 
The signal to noise ratio of the light scattered from a 
spherical water droplet of known size is used as a measure 
60 of the ratio of light scattered from the glare spot to light 
scattered by the other part of the sphere surface. This 
information is then used to predict the shape and size of the 
reflection glare spot. Hence, information contained only in 
the reflection glare spot has been related to the droplet size. 
65 However, until now particle sizes (diameters) have not been 
determined as a function of the distance between a reflection 
glare spot and a transmission glare spot. 
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The present invention provides a new and improved FIG. 6 illustrates glare spot intensity fields for a range of 

apparatus and method which overcomes the above- particle diameters, 

referenced problems and others. 


SUMMARY OF THE INVENTION 

A method for determining a mass flux of an entrained 
phase in a planar two-phase flow records images of particles 
in the two-phase flow. Respective sizes of the particles are 
determined as a function of a separation between spots 
identified on the particle images. Respective velocities of the 
particles are determined. The mass flux of the entrained 
phase is determined as a function of the size and velocity of 
the particles. 

In accordance with one aspect of the invention, the 
recording step includes recording an image of a transparent 
particle. 

In accordance with another aspect of the invention, glare 
spots are identified on the particle. The particle size is 
determined as a function of a separation between the glare 
spots. 

In accordance with another aspect of the invention, the 
step of determining the velocity includes determining the 
velocity as a function of a velocimetry of the particles within 
the images. 

In accordance with a more limited aspect of the invention, 
the step of determining the velocity as a function of the 
velocimetry includes obtaining two exposures of the respec- 
tive glare spots and measuring a displacement between the 
two exposures during a specified time interval. Either the 
average displacement (correlation processing) or individual 
particle displacements can be estimated depending on the 
data processing technique applied and the concentration of 
seed particles. 

In accordance with another aspect of the invention, the 
step of determining the velocity as a function of the veloci- 
metry includes detecting a Doppler shift of light. 

One advantage of the present invention is that it makes in 
situ measurements without disturbing a flow field. 

Another advantage of the present invention is that it 
provides a non-intrusive optical technique that offers instan- 
taneous measurements of particle size and velocity at mul- 
tiple spatial points in planar (2-D) fields, thereby permitting 
an estimate of mass flow to be obtained. 

Still further advantages of the present invention will 
become apparent to those of ordinary skill in the art upon 
reading and understanding the following detailed descrip- 
tion of the preferred embodiments. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention may take form in various components and 
arrangements of components, and in various steps and 
arrangements of steps. The drawings are only for purposes 
of illustrating a preferred embodiment and are not to be 
construed as limiting the invention. 

FIG. 1 illustrates scattering angles, 0 and O, used in a Mie 
Scattering algorithm; 

FIG. 2 illustrates a scattered intensity as a function of 0 
for a 100 /<m diameter PSL sphere in water with 532 nm 
wavelength illumination; 

FIG. 3 illustrates the production of transmission and 
reflection glare spots; 

FIG. 4 illustrates glare spot separation as a function of 
particle diameter; 65 

FIG. 5 illustrates an algorithm for scanning PIV images 
and determining the size of detected particles; and 


DETAILED DESCRIPTION OF THE 
PREFERRED EMBODIMENTS 

A technique making use of the characteristics of the glare 
spots predicted by Mie theory is used in the present inven- 
tion for particle sizing using Digital Particle Image Veloci- 
metry (DPIV). A DPIV system is used to image the light 
scattered from spherical particles under controlled condi- 
tions at an observation angle of 90°. This recorded signal is 
representative of the Fourier transform of the Mie scattered 
electric fields incident on the lens of the imaging system. An 
evaluation of the present invention has been performed via 
a computer simulation program, which incorporates Mie 
theory and the Fourier transform properties of lenses on the 
scattered electric field. The location of the glare spots 
resulting from reflected and transmitted rays in the direction 
of observation are an indicator of particle size. While this 
measurement technique is not as accurate as point measure- 
ment techniques such as Phase Doppler Particle Analysis 
(PDPA), it has the potential for providing reasonable accu- 
racy for size and high accuracy for velocity estimates in 
planar 2-D fields so that mass flow rates can be determined. 
The determination of mass flow rates are essential in many 
production type facilities (e.g., coal furnaces). Experiments 
have been performed to verify the simulation and to deter- 
mine the feasibility of the present technique for measuring 
particle size from Digital Particle Image Velocimetry 
(DPIV) imagery. 

Mie Theory Background 

Mie theory provides the exact solution of electromagnetic 
waves scattered from a spherical particle that includes all of 
the diffraction, specular reflection, interference, and refrac- 
tion effects. At an observation angle of 90°, two (2) glare 
spots are distinguishable on the surface of a transparent 
absorbing particle. The positions of these glare spots are 
described in Mie theory as a function of the particle 
diameter, the relative refractive index, and the angle from 
the optical axis at which the particle is being viewed. 
Knowledge of these positions is used to characterize the 
particles. 

The Mie scattering problem defined here assumes plane 

wave illumination with the electric field E polarized in the 
x-direction (perpendicular to the scattering (y-z) plane) 
incident on a homogeneous, absorbing sphere in a non- 
absorbing medium. This polarization is referred to as 
s-polarization or transverse electric (TE) polarization. With 
reference to FIG. 1 , a particle 10 within a two-phase flow is 
assumed to be in focus at the CCD detector plane. Angles 0 
and $ are defined from the x and z axes to a scattering plane 
12 in which a detector (e.g., CCD camera) 14 is located. E inc , 
B,„ c , and k inc are the electric and magnetic fields and the 
propagation direction of the incident illumination. The 
detector 14 communicates with a processing device 16 , 
which, as discussed below, calculates the size and velocity 
of the particle 10. 

The Mie scattering solution has been derived for the 
problem of the scattering of plane electromagnetic waves 
from a sphere. The magnitude of the scattered intensity is 
defined in Equation 1 as: 
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( 1 ) 


where E 0 is the magnitude of the incident electric field, S 1; 
and S 2 are the scattering amplitudes for the cf> and 0 polar- 
izations of the scattered light, respectively, fi 0 is the perme- 
ability of free space, c is the speed of light, r is the distance 
from the particle to the detector plane, and k is the propa- 
gation constant. This relation shows the polarization depen- 
dence of the scattered electric field from a particle at the 
observation angles 0 and 4>. In the case of PIV imaging, the 
0 range is about 90° and the 0 and <|) range is about 0° or 90°, 
depending on whether the incident electric field is parallel 
(p-polarized, or transverse magnetic (TM) polarized incident 
light) or perpendicular (s-polarized or TE polarized incident 
light) to the scattering plane. When <(>=0°, the Sj contribution 
is negligible since sin 0° is equal to zero (0). Conversely, 
when 4>=90° the S 2 contribution is negligible. For 
s-polarized light (<(>=90°), the reflection and transmission 
glare spots are approximately equal in amplitude and this is 
the preferred illumination state. For p-polarized light (<p=0°), 
the transmission glare spot is significantly dimmer than the 
reflection glare spot, which is not the preferred condition for 
obtaining accurate estimates of the glare spot peak locations. 
FIG. 2 shows how the Sj and S 2 electric fields 20 , 22 , 
respectively, vary over a typical range of 0 values. 

When the incident light is unpolarized, or randomly 
polarized, the magnitude of the scattered intensity is defined 
in Equation 2 as: 

El 1 rl ,1 ,1 (2) 

'^ = 2^cM2' Sm+ 2' Sm \ 


which has equal contributions of S 1; and S 2 scattered electric 
fields. 

The light scattered by the particles is imaged onto a 
detector in order to make a measurement. The work done on 
the electric field imaged onto a CCD array by a lens is 
analogous to taking the Fourier transform of the input 
electric field. Taking the Fourier transform of the electric 
field signal and multiplying it by its complex conjugate 
yields the light intensity distribution on the CCD array (see 
Equation 3). 




scattered) 


(3) 


where is the Fourier Transform, * represents the complex 
conjugate, l scatrere j is the scattered intensity imaged on the 
CCD detector and E scattere j the scattered electric field at the 
lens plane of the optical system. This form of the signal 
makes the scattering effects, such as reflection and 
transmission, visible in the form of glare spot structures in 
the particle image. 

Glare Spots 

With reference to FIG. 3 , glare spots 30 , 32 are bright 
spots resulting from off-axis light scattering by a spherical 
particle 34 . According to ray theory a glare spot is seen at the 
exit points of any rays reflected or refracted in the direction 
of the observer. FIG. 3 demonstrates ray theory for the 
reflected rays 36 and transmitted rays 38 at a 90° observation 
angle. The points 42 and 44 correspond to the reflected and 
transmitted glare spots 30 , 32 , respectively. The ray denoted 
as the reflection ray 36 hits the surface of the spherical 
particle 34 and is reflected from the surface in the direction 
of observation. The ray denoted as the transmission ray 38 
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is internally reflected once and is then transmitted out of the 
particle 34 in the direction of observation. 

The distance Xq+Xj between the glare spots 30 , 32 is 
directly related to the diameter of the particle 34 . Therefore, 
5 it is possible to use the glare spot characteristics to determine 
particle/droplet sizes. As the particle diameter increases so 
does the separation distance between the glare spots. The 
glare spots 30, 32 are oriented in the direction of propagation 
of the incident light. The location of glare spots is calculated 
10 based on the size parameter (2 ita/X), where a is the particle 
radius and X is the illumination wavelength), and the obser- 
vation angle, 0 O , by studying the Fourier transform of the 
geometric scattering amplitudes and the Mie scattered elec- 
tric fields. The essence of the invention consists of estimat- 
15 ing the particle size from a particle image on a CCD detector. 

The glare spot positions relative to the particle centroid, 
x„„ in the object plane are determined and then related to the 
size parameter and observation angle. Therefore, if the 
positions are known at a specific observation angle, the size 
20 of the particle may be determined. The positions of the 
reflection and transmission glare spots on the detector plane 
are defined as: 

.. % (4) 

x a - -aM cos — 

25 2 

9 ^ ‘ (5) 

A‘i = na Msin — + 1 - 2«cos -2- J 

30 where n is the ratio of the index of refraction of the particle 
material to the index of refraction of the medium, a is the 
particle radius, and M is the optical system magnification. 
The measured separation between the glare spots on the 
CCD detector leads to determination of the particle diameter 
35 A p via the following formula: 

. 2A£p (6) 



. B 0 

1 9 °\ 

| — m COS — | + 

nm sin— - 
2 

n 2 + 1 -2ncosy 


where A is the number of pixels separating the glare spots, 
45 f. p is the CCD pixel size, M is the optical system 
magnification, n is the relative index of refraction, and 0 O is 
the observation angle. This relation is used to analyze data 
obtained from PIV images to estimate particle sizes in a flow 
field. 

50 Using Glare Spot Separation to Estimate Particle Diameter 
A one-to-one correspondence exists between the particle 
size and the glare spot separation according to Mie theory. 
Equation (6) gives the particle diameter in micrometers as a 
function of glare spot separation. A, in pixels, CCD pixel 
55 size, e /? , in micrometers per pixel, relative index of 
refraction, n, optical system magnification, M, and obser- 
vation angle, 0 O . FIG. 4 shows a plot 50 of glare spot 
separation as a function of particle diameter based on this 
equation for a polystyrene latex sphere in water ranging in 
60 diameter from 45 fim to 160 /«n with M=0.31, e. p =6.7 
//m/pixel, n=1.195, and 0 o =9O°. The separation with respect 
to particle diameter is linear showing a one-to-one relation- 
ship between the glare spot separation and the particle 
diameter. Therefore, the glare spot separation is a viable 
65 indicator of the particle size. However, due to the low spatial 
resolution of a CCD camera, it is still difficult to extract the 
particle sizes to high accuracy based on the spot separation. 
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The accuracy of the particle size estimates is determined 
by the accuracy by which the two (2) glare spot peaks are 
determined. It has been shown that the uncertainty in esti- 
mating the center location of a Gaussian signal is propor- 
tional to the ratio of the width of the signal to the square root 
of its amplitude. In standard PIV data processing, the 
particle images or correlation peaks are typically estimated 
within an accuracy of about 0.2 pixels, provided the particle 
images span 1-2 pixels across their width. The glare spots 
from the particle images satisfied this criterion. Hence, the 
glare spot separation is estimated to an accuracy of about ^2 
(0.2), providing particle size estimates with accuracies on 
the order of ^2 (0.2)/A. The error in estimating particle size 
is inversely proportional to the number of pixels separating 
the glare spots. Therefore, the CCD camera used for such 
measurements must have sufficient resolution (high pixel 
count) to distinctly image the glare spots. 

Description of Planar Size and Velocity Measurement Tech- 
nique 

A technique making use of the characteristics of the glare 
spots predicted by Mie theory is used in the present appli- 
cation of simultaneous measurement of particle size and 
velocity using PIV. A PIV system is used to image the light 
scattered from spherical particles at an observation angle of 
90°. This recorded signal is representative of the Fourier 
transform of the Mie scattered electric fields incident on the 
lens of the imaging system. As shown in the previous 
section, the location of the glare spots resulting from 
reflected and transmitted rays in the direction of observation 
is a good indicator of particle size. Recall that it is crucial to 
have sufficient spatial resolution to image the glare spots. 
This becomes very important with small particles where the 
glare spot separation is minimal. The centroids of the glare 
spot peaks, and hence the separation between the glare spots, 
are determined using a three-point Gaussian peak estimate 
technique. Using Equation (6), the particle diameter is 
calculated based on the measured glare spot separation and 
the optical system parameters. Software for performing the 
size estimations is described in the next section. Velocity 
estimates are obtained via PIV correlation processing per- 
formed on the same images using any commercial or custom 
PIV processing software. Once the size and velocity are 
obtained, the mass flux of the particle is determined accord- 
ing to conventional techniques. Alternatively, particle track- 
ing velocimetry techniques can be employed, wherein the 
individual particle velocities are determined, which can be 
directly matched to the individual particle size estimates 
obtained from the glare spot separation. 

Sizing 

With reference to FIG. 5, an algorithm 60 has been 
developed to scan the PIV images for particle images and 
determine the size of each particle detected. An image file is 
read in a step 62. Thresholding is applied within the image 
(e.g., global thresholding), in a step 64, to reduce noise 
within the image. Global thresholding imposes a user- 
defined value that limits the minimum non-zero intensity 
value of the image. Any pixels below the threshold value are 
set to zero (0). The thresholding eliminates background 
noise in the image, thereby making it easier to discriminate 
particle images from the background and aiding in the 
separation of overlapping particles. It is important not to 
pick a threshold level so high that the two (2) glare spots are 
separated from each other. Alternatively, a local adaptive 
thresholding technique could also be employed. Local 
thresholding would offer better image extraction than a 
global technique, since variations in background levels in 
the image typically require a higher than optimum global 
threshold level setting. 
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The image pixel array is scanned in a step 66. Preferably, 
the image is scanned from left to right, and top to bottom for 
non-zero (0) pixels. Whenever a non-zero pixel is found, the 
eight (8) adjacent pixels are scanned and any non- zero pixels 
5 are grouped along with the initial pixel to form a particle 
image array. Then the eight (8) adjacent pixels to each of the 
non-zero pixels are scanned and non-zero pixels are grouped 
into the same particle image array. Pixels cannot belong to 
more than one (1) of the particle image arrays. 

10 Scanning continues until no more adjacent non-zero pix- 
els are found and the next particle image is found. Each 
particle image array is analyzed to determine if the image is 
saturated (i.e., more than one (1) pixel in the array has the 
maximum intensity value as determined by the dynamic 
IS range of the camera). The sizing algorithm does not accom- 
modate saturated images since it is difficult to accurately 
determine the locations of the glare spots within the particle 
image when they are saturated. Therefore, saturated particle 
images are rejected. 

20 The glare spot positions and particle diameter are deter- 
mined in a step 68. If the particle is not saturated, the 
maximum intensity in the image is found. The values of the 
pixels above and below, and to the right and left of the 
maximum intensity pixel are used to calculate the 3-point 
25 Gaussian estimate of the glare spot centroid. Next, the 
particle image array is scanned to the left and right (or top 
and bottom in the case of a vertical propagating light sheet) 
to find the point where the slope of the intensity distribution 
curve changes sign, indicating the start of the second glare 
so spot. The particle image array is then scanned from the point 
of slope change to the edge of the array, continuing in the 
same scanning direction, in search of the maximum intensity 
within that area of the array. This maximum intensity 
indicates the location of the second glare spot. The same 
35 3-point Gaussian estimate approach is taken to estimate the 
second glare spot centroid. The difference between the 
centroid positions represents the variable A in Equation 6 for 
calculating the particle image diameter. This same process is 
followed for every particle image array within the entire PIV 
40 image. 

If no slope change is found, then a second glare spot does 
not exist and that particle is rejected. The particle may be too 
small to distinguish between the two (2) separate glare spots 
and, therefore, sizing is not possible using this method. The 
45 locations of the sized particle images and the calculated 
diameter are output to a file. 

Particle diagnostics (e.g., diameter, position, distribution 
mean, and standard deviation) are output to a processor in a 
step 69. 

50 FIG. 6 illustrates glare spot intensity fields for a range of 
particle diameters. Note that intensity fields 70a, 70ft, 72a, 
72 ft, 74a, 74ft, 76a, 76 ft, 80a, 80ft, 82 are very well defined 
for particles having relatively larger diameters. For example, 
the glare spot intensity fields 70a, 70ft for a particle having 
55 a diameter of 150 microns are much more defined than the 
intensity fields 80a, 80ft for a particle having a diameter of 
1 micron. The intensity field 82 only indicates a single glare 
spot and, therefore, would be rejected by the algorithm 
described above. 

60 The preferred embodiment of the present invention has 
been described with respect to estimating the particle size as 
a function of glare spot separations. However, it is also 
contemplated to determine particle sizes as a function of a 
ratio of the amplitude of the glare spot intensities. 

65 Furthermore, it is also contemplated to estimate the diam- 
eter of the particle by iteratively computing the electric field 
scattered by a particle and determining the resultant electric 
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field detected by the CCD sensor. The process is repeated 
until the correct particle size is determined. Convergence to 
the correct particle size is determined when the computed 
image of the electric field on the CCD sensor matches the 
measured glare spot separation and intensity distribution for 5 
a given particle image. 

The invention has been described with reference to the 
preferred embodiment. Obviously, modifications and alter- 
ations will occur to others upon reading and understanding 
the preceding detailed description. It is intended that the 10 
invention be construed as including all such modifications 
and alterations insofar as they come within the scope of the 
appended claims or the equivalents thereof. 

Having thus described the preferred embodiment, the 
invention is now claimed to be: 15 

1. A method for determining a mass flux of a particle, 
comprising: 

recording an image of the particle entrained in a two- 
phase flow, using a camera; and 
using a processor for determining a size of the particle as 
a function of a separation between spots identified on 
the particle, determining a velocity of the particle, and 
determining the mass flux of the particle as a function 
of the size and velocity, wherein the spots are glare 25 
spots and the separation between the glare spots is 
determined as: 
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where A p is an estimate of the particle diameter, n is a ratio 
of an index of refraction of a material of the particle to an 
index of refraction of a medium, a is a radius of the particle, 

M is an optical system magnification, A is a number of pixels 
separating the glare spots on a surface of a CCD, e p is a size 45 
of the pixels in the CCD, and Q a is an observation angle. 

2. The method for determining a mass flux of a particle as 
set forth in claim 1, wherein the recording step includes: 

recording an image of a transparent particle. 

3. The method for determining a mass flux of a particle as 50 
set forth in claim 1, further including: 

identifying glare spots on the particle, the particle size 
being determined as a function of a separation between 
the glare spots. 

4. The method for determining a mass flux of a particle as 55 
set forth in claim 1, wherein the step of determining the 
velocity includes: 

determining the velocity as a function of a velocimetry of 
the particles within the images. 

5. The method for determining a mass flux of a particle as 
set forth in claim 4, wherein the step of determining the 
velocity as a function of the velocimetry includes: 

obtaining two exposures of the respective glare spots of 
the particles entrained in the fluid; and 6S 

measuring a displacement between the two exposures 
during a specified time interval. 
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6. The method for determining a mass flux of a particle as 
set forth in claim 4, wherein the step of determining the 
velocity as a function of the velocimetry includes: 

detecting a Doppler shift of light. 

7. An optical flow meter for determining a mass flux of a 
particle, comprising: 

a camera for recording an image of the particle entrained 
in a two-phase flow; and 

a processor for determining a size of the particle as a 
function of a separation between spots identified on the 
particle, determining a velocity of the particle, and 
determining the mass flux of the particle as a function 
of the size and velocity, wherein the spots are glare 
spots and the separation between the glare spots is 
determined as: 
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where d p is an estimate of the particle diameter, n is a ratio 
of an index of refraction of a material of the particle to an 
index of refraction of a medium, a is a radius of the particle, 
M is an optical system magnification, A is a number of pixels 
separating the glare spots on a surface of a CCD, e^, is a size 
of the pixels in the CCD, and Q a is an observation angle. 

8. The optical flow meter for determining a mass flux of 
a particle as set forth in claim 7, wherein a Gaussian peak 
location estimate is used for determining a location of 
respective peaks of the glare spots, the separation between 
the glare spots being determined as a function of the 
locations of the peaks. 

9. The optical flow meter for determining a mass flux of 
a particle as set forth in claim 7, wherein the camera is a 
CCD camera. 

10. The optical flow meter for determining a mass flux of 
a particle as set forth in claim 7, wherein the particles are 
transparent. 

11. A method for determining a size of a particle, the 
method comprising: 

receiving an image of the particle entrained in a two- 
phase flow into a processor; 

reducing background noise within the image; 

grouping the pixels having non-zero values into respec- 
tive particle image arrays; 

identifying glare spots within the image as a function of 
the particle image arrays; and 

determining the size of the particle as a function of a 
separation between the glare spots, wherein the sepa- 
ration between the glare spots is determined as: 
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where d is an estimate of the particle diameter, n is a ratio 
of an index of refraction of a material of the particle to an 
index of refraction of a medium, a is a radius of the particle, 

M is an optical system magnification, A is a number of pixels 
separating the glare spots on a surface of a CCD, e p is a size 15 
of the pixels in the CCD, and 0 O is an observation angle. 

12 . The method for determining a size and a velocity of 
a particle as set forth in claim 11, wherein the reducing step 
includes: 

limiting non-zero intensity values of pixels within the 20 
image. 

13 . The method for determining a size and a velocity of 
a particle as set forth in claim 12, wherein the limiting step 
includes: 

determining a global threshold intensity value for the 25 
pixels within the image; and 

setting intensity values of pixels below the global thresh- 
old to zero. 
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14 . The method for determining a size and a velocity of 
a particle as set forth in claim 13 , further including: 

determining a local threshold for discriminating the par- 
ticle within the image. 

15 . The method for determining a size and a velocity of 
a particle as set forth in claim 11, wherein the grouping step 
includes: 

scanning the image for the pixels having the non-zero 
values; 

identifying one of the pixels as having the non-zero value; 

identifying pixels adjacent to the pixel having the non- 
zero value; 

grouping any of the adjacent pixels having the non-zero 
values into the particle image array; 

identifying subsequent pixels adjacent to each of the 
adjacent pixels having the non-zero value; and 

grouping any of the subsequent pixels into the particle 
image array. 

16 . The method for determining a size and a velocity of 
a particle as set forth in claim 11, further including: 

rejecting ones of the particle image arrays that are satu- 
rated. 



